OKHEP-94-13 

Absence of Species Doubling in Finite-Element Quantum 

Electrodynamics 



Kimball A. Milton* 

Department of Physics and Astronomy, University of Oklahoma, Norman, OK 73019 USA 
G\ . (February 1, 2008) 

o 

Q. 

q\ ; Abstract 

i ! 

In this letter it will be demonstrated explicitly that the finite-element formu- 

O . 

£Sj . lation of quantum electrodynamics is free from fermion doubling. We do this 

m _ 

by (1) examining the lattice fermion propagator and using it to compute the 
. one-loop vacuum polarization on the lattice, and (2) by an explict computa- 
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tion of vector and axial- vector current anomalies for an arbitrary rectangular 



Qh. lattice in the Schwinger model. There it is shown that requiring that the 



vector current be conserved necessitates the use of a square lattice, in which 
case the axial- vector current is anomalous. 

I. INTRODUCTION 

The application of the finite-element method to the Heisenberg equations of motion of 
a quantum theory was proposed over a decade ago 0. It was immediately recognized that 
the formulation was unitary, in that canonical commutation relations are maintained at 
each lattice site. (This was in addition to the fact that of all possible discretizations of the 
equations of motion, the finite-element prescription is the most accurate 0.) The Dirac 
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equation was studied shortly thereafter , and it was found that the dispersion relation did 
not admit fermion doubling, usually present in lattice formulations. The reason that the 
present formulation could evade the no-go theorems [|J is that no local Hamiltonian in the 
Schrodinger sense exists. However, this conclusion has remained somewhat controversial, 
and it is the purpose of the present letter to offer new evidence for the absence of species 
doublers. In particular in Sec. [ij we display the free fermion propagator, and show that only 
a single fermion is represented, and that when these propagators are used in a dynamical 
loop calculation, the anomaly in the Schwinger model is indeed recovered. In Sec. |(| we 
explicitly compute the divergence of the vector current in (1 + 1) dimensions, and show 
for a square lattice that the vector current is conserved, while the axial-vector current is 
anomalous. We repeat the calculation for small rectangular lattices, and show that only for 
square lattice case can a vector anomaly be avoided. 

Abelian and non-Abelian gauge theories in the finite-element context were discussed, 
respectively, in || and |J. The Schwinger model was treated in this context previously in 
|5|,|7]]. A review of the entire program appears in ||. 

II. FREE LATTICE FERMION PROPAGATOR. ONE-LOOP CALCULATION OF 

VACUUM POLARIZATION 

From the free finite-element lattice Dirac equation, 



where /i is the electron mass, h is the temporal lattice spacing, A is the spatial lattice 
spacing, m represents a spatial lattice coordinate, n a temporal coordinate, and overbars 
signify forward averaging: 
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(2.1) 




(2.2) 



it is easy to derive the free fermion propagator ||, 
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p 

7° sin ^ + (/i — 7 • p)/i cos 2 
cos h(Q — ze) — cos hVt 

Here L = MA is the length of the spatial lattice, and 

2t 

p = — , uj = p° = Vp 2 + ^ 2 , (t p )i = tan^7r/M. (2.4) 

The mass-shell "energy" is defined in terms of A, the eigenvalue of the Dirac transfer 
matrix, 

A = 1 + = e «W\ (2.5) 

1 - ihu/2 v ; 



Notice that we may solve ( |2.5| ) for u: 



2 Ml , s 

lo = — tan — . (2.6) 
h 2 

We have taken M, the number of lattice points in a given spatial direction, to be odd, so 
that ij} is periodic on the spatial lattice. 

We note that the singularities in the fermion propagator fl2.3p are unique: That is, for a 
given lattice momentum p there is a pole at a single energy Cl between and ir/h. A species 
doubler would be signalled by having the minimum value of u, /i, occur at p = 0(mod M/2) 
not 0(mod M). 

We confirm this conclusion by presenting the results of a dynamical look calculation 
based on the propagator ( |2.3| ) of the photon polarization tensor in the Schwinger model 
(/i = 0, d = 2) [For details of this calculation, see 0.] In the continuum, the polarization 
tensor is transverse, 

/ h^h u \ 

n^ = H(A; 2 ) , (2.7) 
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where an easy calculation shows that 11(0) = — e 2 /7r, which gives the boson mass in the 
model. It is straightforward to carry out the corresponding calculation using the lattice 
propagator; the result is shown in Figure 1. The Schwinger model corresponds to fi = 0. 



In fact, it is clear that for very small v = //A the lattice calculation diverges, which is as 
expected when v ~ M -1 . There is, for all values of r = h/A, a strong shoulder at the 
continuum value of 11(0). (The sensitivity as r — ► 1 is a lattice artifact.) Such behavior is 
completely consistent with the absence of species doubling. 

III. INTERACTIONS. AXIAL- VECTOR ANOMALY 

Interactions of an electron with a background electromagnetic field is given in terms of 
a transfer matrix T: 



if) n+1 = T n ijj n , 



(3.1) 



which is to be understood as a matrix equation in m. Explicitly, in the gauge A = 0, |T0 



where 



T = 2U- 1 - 1, U = 1 + 



ih/j.'y h 
~ A 



7°7 • 2>, 



(3.2) 
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(3.3) 



Here 

1, m > m', 
0, m = m', 
—1, m < m', 

and (the following are local and unaveraged in m^, n) 
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(3.4) 



(3.5a) 



and 
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sgn (mj - m,)sgn (m" - mJ)Cn»y> 



mv=l 



(3.5b) 



with 



sgn (m - rri) = e m>m > - b m)m <- (3.6) 

Because T> is anti-Hermitian, it follows that T is unitary, that is, that <p m ,n = ipm,n is the 
canonical field variable satisfying the canonical anticommutation relations 

{0m,n>0L>} = ^m,m'- ( 3 - 7 ) 

In this letter we will consider the Schwinger model, that is, the case with dimension 

d = 2 and mass p, = 0. Because the light-cone aligns with the lattice in that case, we first 

set h = A. Then we see that the transfer matrix for positive or negative chirality, that is, 

for the eigenvalue of 275 = 7°7 1 equal to ±1, is 

1±V , . 

T± = . 3.8 

1=FX> v ' 

From (|3.3| ) we see that the numerator of T is 

(1 ± V) m , m , = [5 m ^e^ =F (-l) m+m 'e m , m ,e--"^—'] sec C, (3.9) 
while it is a simple calculation to verify that the inverse of this operator is 

(! + V)m,mf = g i 6m ' m ' + ^, m '-ie~ 2 ' Cm ') , (3.10a) 

(1 " ^)mU' = 2 ^ m ' m ' + ^'+l e2<m ) • ( 3 ' 10b ) 

It is therefore immediate to find 

(T+) m , m / = 5 m , m / +1 e 2iCm , (3.11a) 

(^-)m,m' = ^m+l, m 'e~ 2<m ', (3.11b) 

which simply says that the + (— ) chirality fermions move on the light-cone to the right (left), 
acquiring a phase proportional to the vector potential. Solution (|3.11|) directly implies 
the chiral anomaly in the Schwinger model, as shown in |J. We will now expand on and 
generalize that calculation. 

First, the canonical anticommutation relations (|3.7|) show that we can write the mo- 
mentum expansion of the (free) Dirac field as (if M is even we replace p by p + 1/2 in the 
exponent) 



s.p 



£ ^(b p , s u p , s \- n e* p - m2 ^ M + dl s v p , s X n e- tpm2 ^ M ) } (3.12) 



where the spinors are normalized according to 

J2u ± uW = = ±A ± , (3.13) 

s 

where u = «_ and v = u+. The creation and annihilation operators satisfy (d = 4) 

{ b P,s' b l',s'} = ^W^,^ (3.14a) 
{ d p^ d t',s'} = ^W^'i (3.14b) 

and all other anticommutators of these operators vanish. It is then easy to verify that for 
M even, for the particular case of d = 2 and \i = 0, that in the Fock-space vacuum^ 

^& = T TL ^-_i >/M , (3.15a) 



while for M odd, 



,(±)U(±) v = j cos(m - mQTr/M - (-1)' 
m ' n ?m ' + 2L sin(m-m>/M 



(3.15b) 



In both cases, the vacuum expectation value is taken to zero if m = m' . 

It is then quite straightforward to compute the vacuum expectation value of the finite- 
element divergence of the the vector current. On the lattice the gauge-invariant current is 
written in terms of the all- averaged field 

^m,n = ^m,n = ^(0m,n+l + 4>m,n), (3.16) 

rather than the canonical field <fi. That is, the current is 

JL, n = e*L,„7V* m> n = + Ij)]»7°y [(1 + T n )<P n ] m . (3.17) 



l In the following we omit terms arising from reordering the operators on the left-hand side of 
( 3.15| ) using the anticommutator ( |3.7| ); these are local in m, m', and hence do not affect the r = 1 



results. In general they only contribute an extraneous spatial constant to (J°). 



(For a discussion of why this choice of current is used, see the Appendix of ||.) In the 
absence of interactions it is easy to show that 

( u dpJ fi ») = 0, (3.18) 

where the quotation marks signify a finite-element lattice divergence and the brackets rep- 
resent a Fock-space vacuum expectation value. In two dimensions, that divergence is 



+ q * fe+l,n+l'75fe+l,?i+l + fe+l,n'75fe+l,n '4 , m,n+l''l^m,n+\ ^Tn^l^mjl) ■ (3.19) 
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For the case of a square lattice, h = A, this is simply expressed in terms of eigenvectors of 
ij5, which in turn may be expressed in terms of the canonical field 4> m ,n+i at the intermediate 
time. A short calculation reveals that 

= + (3.20) 

where 



a(±) - JL 

A ~ Ah 



(,(±)t A±) „-i(Cm+l,w+Cm+l,w+l) 
; m,n+lr m +l,n+l c 



,(±)f ,(±) i(C ra +i,„ + Cm + l,n+l) 

; m+l,n+lVm,n+l c 



2i sin((^ jn+lin — £ m+lira+1 

) (3.21) 



The vacuum expectation value of cancels that of A 1 - \ and hence the vector current is 
conserved, because (q is the eigenvalue of ry 5 ) 

Mi ^ « cos2 7r / 2M 

(tfirf) = I ^ Usctt/M, (3.22) 

for M odd or even respectively. On the other hand, the divergence of the axial-vector 
current, 

"«9 M Jf = A {+) - A { ~\ (3.23) 

is not zero: 



<"<Vf) 



2e 1 
hL sin7r/M 



cos 2 7r/2M 



sin(C m+ i jra+ i — Cm+l,n) COs(Cm+l,n+l + Cm+l,n)- 



(3.24) 



In fact, expanding this in powers of ehA, we find (because E — A) 



[VP) 



cos 2 tt/2M 



M sin 71 /M 



E (1 + +0{{ehA) 2 )) 



(3.25) 



For M = 2 the error in the coefficient relative to the continuum value of the coefficient, e 2 /7r 
is about 50%, while for M = 3 the error drops to about 10%; in general, the relative error 
is of order M~ 2 . 

We now consider a rectangular lattice, h ^ A. Because this case is rather more compli- 
cated that the square lattice considered above, we content ourselves with the two smallest 
possible lattices, M = 2 and M = 3. For the M = 2 case the transfer matrix 

l±rV 



(3.26) 



where r = h/A, is easily found to be (the subscript on ( indicates the spatial coordinate) 



/ 



T- 



± 



(l-r 2 )(l + e 2 <) T 4 



re 



2iCi 



\ 



where 



^ ±4re 2 ^ (l-r 2 )(l + e 2 <) y 



D± = e 2K {lTr) 2 + (l±rf 



(3.27) 



(3.28) 



Note that 



D = D* ± D ± = 2[1 + 6r 2 + r 4 + (1 - r 2 ) 2 cos 2(]. 



(3.29) 



A straightforward calculation now reveals (the superscript now indicates the time coordinate) 
that 

4er(sin2Cj 2) -sin2Cf } ) 



(A) 



-<4 2 > = 



(3.30) 



which are expressed in terms of potentials at the intermediate time 2, while 

(J°n) = -(J°2i) = -(J°uf 1 \ (3-31) 



where the last notation means that the expression is the same as (|3.30|) except that the 
potential is at time 1. In the same way, we find 

/ t1 \ / t1 \ 4er 2 (sin2Cf 2) +sin2Cf ) ) 

Wl2/ — W22/ — AD^ [o.OZ.) 

and 

(Jn) = (4> = (^ 2 ) C(1) - (3-33) 
It follows immediately that the vector current is conserved: 

("fyj"") = 0. (3.34) 

However, the axial- vector current is anomalous 

8eh( sin2d 2) sin 2^ A . , 

[ <V 5 ) - ^g" I D{ 2) D (l) J ■ ^- 6b ) 

Note that this reduces to the previous r = 1 result, G3.24| ). 

In the same way the calculation at M = 3 can be carried out. Here the transfer matrix 
is given by 

T± = x (3.36) 

f (1 -r 2 )(l ±r + (1 Tr)e 2i< ) 4r(r T l)e 2i(Cl+Cs) 4r(r ± l)e 2iCl ^ 

4r(r±l)e 2 ^ 2 (1 - r 2 )(l ± r + (1 =f r)e 2iC ) 4r(r =F l)e 2i( ^ 1+ ^ 2 ) 

I 4 r ( r T i) e 2i(C 2 +C3) 4r(r±l)e 2i?3 (1 - r 2 )(l ± r + (1 =f r)e 2iC ) ^ 

where 

D± = (lir) 3 + (lTr)V ,( . (3.37) 
Here the common denominator is 



D = D* ± D ± = 2[1 + 15r 2 + 15r 4 + r b + (1 - r 2 fcos2C]. 



(3.38) 



We find 



(A) 



4er 



(1 + 3r 2 ) (sin 2^ - sin 2(f) + (1 - r 2 ) (sin2(C 1 (2) + (P) 



-sin2(Cf ) + Cf: 



(3.39) 



where the corresponding expressions for ( J^) and (J| 2 ) are obtained by translation. Again, 
(Jii) is obtained by replacing by 



(J°n) = -(J' 
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(3.40) 



The vacuum expectation value of J 1 is 

4er 2 



(A) 



"(3 + r 2 ) (sin2C 1 (1) + sin 2^) + (1 - r 2 ) (sin2(cj 1) + C^) 



+ sin2(d 1) + d 1) )-2sin2Ci 1) 



and (Ji 2 ) is obtained by replacing by Q 2 \ 



(3.41) 



(3.42) 



Now, in general, neither vector nor axial-vector current is conserved: 
/ tm» \ _ J_ / 70 , 70 _ 70 _ t0 \ , J_ / rl , rl _ t1 _ rl \ 

\ UfiJ / — 2/^^ 22 12 21 n ' 2A ^ 22 21 12 n ' 

= ^§ [(1 - r 2 ) (sin2 Cl (2) - sin2 C f) + (1 + r 2 ) (sin2( Cl (2) + g>) 

-sin2(Cf + d 2) ))l +(C (2) -C (1) ), (3.43) 



> ^^5* ) ~~ 2h^ 22 "^ 12 "^ 21 2A^ 22 ^ 21 ^ 2 "^ n ^ 



4er 



v^Z^A 2 



4(1 + r 2 ) sin2Cf + (1 - r 2 ) (sin2(C{ 2) + (f) + sin2(Cf } + Cf } ) 



-(C (2) -C {1) ). 



(3.44) 



However, if we require the vanishing of the vector anomaly, that is, electric current conser- 
vation, we must set r = 1 and then the result in ( p.24|) follows. 
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IV. CONCLUSION 



In this letter we have examined the axial-vector anomaly in two-dimensional electrody- 
namics in the context of the finite-element lattice. The results seem quite definitive: In 
agreement with arguments based on the dispersion relation (or the fermion propagator), 
chiral symmetry is broken and no fermion doubling occurs. The reader will ask, how can 
this be? The answer lies in the fact that in our Minkowski-space formulation, no Lagrangian 
exists from which the equation of motion (Dirac equation) is derivable. Hence symmetry 
arguments do not imply a corresponding conserved current. This conclusion is supported by 
the results discussed in the Appendix to |J. There, it is shown that in Euclidean space (with 
periodic boundary conditions in all directions) it would be possible to define a Lagrangian, 
from which conserved vector and axial- vector currents could be derived. But, in Minkowski 
space this cannot be done, and the gauge-invariant local currents are anomalous. 
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FIG. 1. Plot of the one-loop vacuum polarization on a rectangular lattice with M = 2533 
spatial lattice sites as a function of v = //A/2. Shown are curves with r = h/A = 0.1, 0.75, 0.99. 
The quantity plotted is -II(0)/(e 2 /7r)- 
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This figure "figl-l.png" is available in "png" format from: 
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